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ABSTRACT 

Context. Broad MgllA2800 and H/3 lines have emerged as the most reliable virial 
estimators of black hole mass in quasars. Which is more reliable? Part of the challenge 
centers on comparing MgllA2800 and H/3 line profiles in order to improve the ± 1 dex 
A/bh uncertainties inherent in single-epoch FWHM measures from noisy spectra. 
Aims. Comparison of MgllA2800 and H/3 profile measures in the same sources and 
especially FWHM measures that provide the virial broadening estimator. 
Methods. Identification of 680 bright Sloan Digital Sky Survey DataRelease 7 quasars 
with spectra showing both MgllA2800 and H/3 lines, at redshift 0.4 < z < 0.75. The 
s/n of these spectra are high enough to allow binning in the "four-dimensional (4D) 
eigenvector 1" optical plane and construction of high s/n composite spectra. 
Results. We confirm that MgllA2800 shows a profile that is ~ 20% narrower as sug- 
gested in some previous studies. FWHM measures for Population B sources (i.e., with 
FWHM of H/3 larger than 4000 km s^^) are uncertain because they show complex pro- 
files with at least two broad-line components involving a nearly unshifted broad and 
redshifted very-broad component. Only the broad component is likely to be a valid 
virial estimator. If H/3 and MgllA2800 are not corrected for the very broad compo- 
nent then black hole mass Mbh values for Population B sources will be systematically 
overestimated by up to A logMeu ~ 0.3 — 0.4 dex. We suggest a simple correction 
that can be applied to the majority of sources. MgllA2800 is the safer virial estimator 
for Population B sources because the centroid shifts with respect to rest frame are 
lower than for H/3. In the broad and very broad component profile interpretation 
this is a consequence of the lower very broad to broad component intensity ratio for 
MgllA2800. Eigenvector-based studies show that effective discrimination of black hole 
mass and Eddington ratio at fixed redshift is not achieved via luminosity binning but 
rather by binning in a "4D eigenvector 1" context that refiects different broad line 
region geometry/kinematics likely driven by Eddington ratio. 

Key words, quasars: general — quasars: emission lines — line: profiles — black hole physics 



1. Introduction 



Estimations of black hole mass (A/bh) Siud Eddington ratio ( L/LEdd) for large and diverse 



samples of quasars (jMarziani fc Sulentic 



2012; 



tio ( j 

ioia . 



Shed |2013| . and references therein) are 
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needed for both astrophysical and cosmological studies. Until recently the majority of 
virial black hole A/bh estimates have not taken profile diversity into account. Rather, 
broad-line FWHM measures have been assumed to be valid virial estimators in all quasars. 



from different lines that often show quite different profile shapes (e.g., 


Barthel et al. 


1990|: 


Baskin & Laor 


2OO5I: Brotherton et al. 


1994 


: Ho et al. 2OI2I: 


Wills et al. 1993: Z 


amfir 
3 and 


et al. 2010). The goal of the present oaoer is to svstematize differences between H/ 



MgiiA2800 in the process of comparing the relative merits of FWHM H/3 and MgiiA2800 
as virial broadening estimators. This is done within the context of a 4D Eigenvector 1 
parameter formalism. 



1.1. The 4D eigenvector 1 parameter space 

Interpretation of quasar spectra in the context of a 4D parameter space stemming from 
eigenv e ctor studies reveals a large d iversity of line profile properties ( Marziani et al. 201ol . 



2003b: Sulentic et al. 



2007 



2000al ) h kely driven by L/I-Eddi with Mbh and source ori- 



entat ion also plaving important roles ( Boroson fc Green ^ 1992 ^ Marziani et al.l 2001 ^ Yip 



et al. I2OO4I ). The 4DE1 (4D Eigenvector 1) formalism owes its inception to a principa l 



component analysis study of systematic trends in quasar spectra (jBoroson fc Greenlll992r ) 
from the the Palomar-Green (PG) survey. The first eigenvector (El) included a correla- 
tion between [Oill] AA4959,5007 and optical Pen strength and an anti-correlation between 



(e-g- . 


Borosonll2002: Grupe 


2004 : Kovacevic et a 


.II2OI0I: Kruczek et al.ll2011: Kuraszkiewicz 


et al. 


2002l:lMarziani et al. 


1996: 


Sulentic et al.l 


2002 


. 2000bl: Tang et all 2012: 


Wang et al. 


2006: 


Yip et al.ll2004). El was eventuallv expanded into 4DE1 parameter space (Sulentic 


et al. 


2OO7I. 2000a) with the addition ClvA1549 profile shift and soft X-rav photon index 



as important parameters that show some of the strongest intercorrelations. The higher 
dimensionality involves parameters that are observatio nally independent an d that corre- 



spond to different phvsical processes (see the reviews bv Sulentic et al.ll2000a and Marziani 
& Sulentic laOli for more detailed discussions). Quasars occupy a well defined sequence in 
the "optical plane" of 4DE1 defined by the PWH M of the H/3 broad component and by 
the i ntensity ratio between PeiiA4570 and H/3 (jBoroson fc GreenI Il992l : ISulentic et al. 
20021 ). The 4DE1 optical plane therefore offers a tool that allows for the definition and 
identification of quasar populations or spectral types. 



1.2. Two main quasar populations 

Exploration of the 4DE1 parameter space gave rise to the concept of two quasar populations 
that, if not truly distinct, are an effective way to distinguish important spectroscopic 
differences between quasars. Pop. A sources show a wide range of broad-line widths 
(FWHM H/3^ 600 - 4000km s~^) with a majority of sources (in low redshift samples) 
between FWHM H/3- 1000 - 4000 km s'^. Few sources with i?Feii> 0.5 show FWHM 
H/3> 4000 km s^^ which was partial motivation for the Pop. A-B formalism. Pop. A 
involves a largely radio-quiet quasar population as few radio-loud quasars (especially with 
double-lobe Fanaro ff-Riley II (FRII) radio morphology) show FWHM H/3<4000 km s"! 
(jZamfir et allboosh . AGN with FWHM H/3< 2000 km s'^ are often referred to as narrow 
line Seyfert 1 (NLSyl) sources but there is continuity in source properties over the full 
Pop. A range of FWHM H/?. The "Population A" definition changes in higher redshift 
samples because the minimum obs erved broad-line FWH M o f H/3 slowly increas es with 
source luminosity (see Figure 11 in Marziani et al. 20091 . and Dultzin et al. 2011 ). Pop. 
A H/3 profiles are best fit with a single symmetric Lorentz function while Pop. B are 
not. Following previous results for H/? we fit all Pop. A MgiiA2800 profiles with Lorentz 
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functions ([Sulentic et alJ l2002l : IVeron-Cettv et al.ll200lh . Pop. B quasars show an even 
wider range of FWHM H /3 than Pop. A sou rces. The full observed range is FWHM H/3 
= 4000 - 40000 km (jWang et allboool ). but values FWHM H/3> 16000 km s"! are 
extremely rare with almost all sources in a low redshift sample between 4000 km s^^ and 
12000 km s~^. Most Pop. B sources show i?FGii< 0.5. Pop B is a mixed radio-quiet and 
radio-loud population with the majority of radio-loud sources within population B. Pop. B 
sources are different from Pop. A because - most fundamentally - they show composite H/3 
profiles. The majority of Pop. B H/? profiles can be modelled with an (almost) unshifted 
broad (BC) and redshifted very broad component (VBC). It is usually assumed that the 
unshifted BC is the same as the single component observed in Pop. A sources while the 
VBC is something non-virial and apparently unique to Pop. B sou rces (both rad i o-quie t 
and radio- loud). Further discussion on Pop. A and B is provided bv lSulentic et al.l (j201lh . 



1.3. Black hole mass estimates 



The equation used to estimate A/bh can be written, under the assumption that line broad- 
ening in quasars is due to virialized motions in the emitting g 



BH 



rBLR(Aw)^ 



G 



(1) 



wher e J- is factor 1 tha t depends on geometry of the emitting gas (e.g.. ICraham et al 



2011 



Onken et al.l 120041 ) , tblr is an emissivity- weighted radial distance of the broad 



line emitting region (BLR), and A?; is the virial broadening term, customarily measured 
from the width of a suitable emission line. The most accurate Mbh estimates co me from 
reverberation mapping measures of rRi R (e.g. lDennev et al I2OO9I J iHorne et al l2004h coupled 



with H/3 rms profile widths (FWHM or vel ocity dispersion cr: 



now have «60 reverberation derived radii (Bentz et al 



201 



Pete rson et al 
2OO9I . 



201 



200J). We 



The cost in 

telescope time for such accurate tblr and rms line width estimates is very high even 
for the brightest sources. It is clear that we must rely on secondary methods in order 
to obtain a large numbers of additional A/bh est imates, as needed for the analysis of a, n 
ever increasing population of catalogued quasars (|Ross et al.ll2012l : ISchneider et ahlboioh . 
These methods involve single-epoch FWHM measures and rBLR estimates derived from the 
apparent correlation between reverberation deri ved radii and source luminosity measures 
(the so-called "Kaspi relation:" iKaspi et al.ll2005l ). We suggest that line measures for large 
source samples can be better utilized if 4D parameter space coordinates for each source 
are taken into account. This is especially true if 4DE1 coordinates are driven by source 
Eddington ratio. 

Important issues connected with the use of Eq. [IJinvolve: 1) the validity of extrapolating 
the Kaspi relation beyond the range of redshift and source luminosity represented in the 
reverberation sample and 2) estimation of the virial broadening Aw. There are two sides to 
virial broadening estimation: a) the choice of a suitable emission line (assumed primarily 
broadened by virial motions) and b) selection of a suitable parameter (FWHM, a) to 
define the virial broadening Aw. Single-epoch FWHM H/3 measures have provided Mbh 
estimates for large samples at low and high redshifts. Measures of the second moment a 
are less useful in the case of a broad lines with more than one component. Beyond redshift 
z ^ 0.7 we lose the most trusted source of single-epoch FWHM measures (i.e. H/3). We 
can either follow H/3 into the infrared and/or adopt FWHM measures of other broad lines 
as virial estimators. We h ave adopted the former approach taking advantage of Very Large 
Telescope (VLT) ISAAC (jMoorwood et al.lllQOsh spectroscopy to follow H/3 through the 
IR windows and obtain accurate as possible Mbh estimates for ? ^50 lu minous quasars in 



the range 2 « 1.0 - 3.0 (jMarziani et al.lboodlsulentic et al]l20od . [20o3 ). We are currently 
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processing a sample of «25 sources where H/3 is redshifted into the K band window (z 
3.0 - 3.7) beyond which H/3 is lost (Sulentic et al. 2013 in preparation). 



1.4. Use ofMguX2800 



The strongest available surrogate broad lines that are redshifted into the optical over a wide 
z range are ClvA1549 and MgiiA2800. The region of Cm] A1909 is a blend of permitted and 
semi- forbidden lines whose apparent intricacy has deterred previous workers although there 



is potential for progress if moderate resolution h igh s/n spectra are available (jMarziani et al 



2011, 


20inl: 


Neerete et al. 


2OI2I: 



(Marziani et al. 19961: Netzer et al. 200?!: iRichards et al. 2011: Shen et al. 20081: 


et al. 2007) 


C1VAI549 also shows a blueshift and profile asymmetry 


in many 


( Brotherton 


1996|: 


Gaskel] 


1982|: 


Richards et al. 


2011; 


Sulentic et al. 


2007 


) challenj 



tic 



validity of the virial assumption. Some consistency can be achieved i f several cauti ons and 
corre ctions are applied to the ClvA1549 measures ( Assef et ahlboil ^ 



20ld and especially in the 4DE1 context as done by 



Denne 



Sulentic et al 



^l20li:G 



reene 



2007h but this 



et al. 

i s not the road to Mb h estimates with uncertainties less than Idex. In the 4DE1 context 
2000bh both H/3 and CivA1549 show strikingly different profile properties 



(Sulentic et al 



for Pop. A and B sources. Such differences cannot be easily quantified with low resolution 
and/or low s/n spectra. 

The rest wavelength of MgllA2800 makes it a potentially valuable ground-based virial 
estimator from z wO.4 up to z w 2.1. Rapid improvements in IR spectroscopy also makes 
MgiiA2800 accessible at higher redshifts. It has b een measured in K b and for the quasar 
with highest known redshift z « 7.085 redshift (jMortlock et al.lboilh . Is MgliA2800 a 
safe H/3 surrogate? Mg+AA2796,2803 is a low ionization line (LIL) doublet long thought 
to be produced in a dense optically-thick r egion where at least part of the Balmer lines 



(jCrandi fc Phillipjll979l : 



Netzer 



19801 ). A possible c omphcation is that MgiiA2800 



is a resonance line and therefore often af fected by broad (jGangulv fc BrothertonI 120081 ) 
and narrow absorption ( Wild et al. 2008h . Several authors have attempted to define a 
relati on between MgllA2800 and H/3 to provide consistent Mbh estimate s from the two 
lines jMcLure fc Jarvisll2002l:lshen fc Liull2012l:lTrakhtenbrot fc Netzeilboid: Vestergaard & 



Peterson l2006riWang et al 
as H/3 because their line widths are correlated. Latest studies however show considerable 
scatter (and deviation from parity) in the FWHM H^ vs. FWHM MgiiA2800 correlation. 
There is strong ev idence that FWHM MgliA2800 is about 20% smaller than FWHM H/3 
(jWang et al.ll2009l ) in many sources. 

We report here on a more detailed comparison of H/3 and MgiiA2800 as virial estimators. 
This work goes beyond several recent works considering line shape properties (shifts and 
asymmetries, presence of a very broad component) that could weaken or invalidate the 
virial assumption. Such a more detailed study requires separate consideration of sources 
with FWHM H/3 les s than or greater th an 4000 km s^^ (Population A and B respectively; 



20091 ). MgllA2800 is assumed to be a virial estimator as reliable 



cf. Pop. 1 and 2 of 



companion paper (jMarziani et al 



Collin et al 



1200^ 



201 



. Results for Pop . A sources were presented in a 
hereafter iPaper if ). The majority (w 80 %) of Pop. 



A sources show an unshifted MgllA2800 profile that is significantly less broad (^20%) than 
H/3 making it the virial estimator of choice. Pop. A H/3 profiles are best fit with a single 
symmetric Lorentz function. Gaussian fits to Pop. A profiles will result in a systematic 
overestimation of Mbh- This paper focuses mainly on Pop. B sources where broader and 
more com plex lines are observed. A subset of Pop. A sources show a significant MgllA2800 
blueshift (jPaper J ) making FWHM H/3 a safer virial estimator. These sources can be 
described as the youngest least-evolved quasars with strong ongoing star formation mong 
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all quasars in the present sample fe.g.. lMatlnii]l20od : [Sani et allboidlSulentic et alJboOO 



Wang et alJbooeh . Extreme values of i?Feii appear to be an empirical signature of quasar 



youth perhaps reflecting extreme low ionization and metal enrichment as well fc.f. Paper 
I). 

The advent of the Sloan Digital Sky Survey (SDSS) makes possible accurate direct 
comparison of H/3 and MgiiA2800 using the brightest 600+ quasars where both lines appear 
in the spectra (^describes the sample of SDSS quasars). Individually most of these spectra 
are too noisy (typically s/n < 5 in the continuum near the two lines) to yield good measures 
and, therefore, a reliable application of A/bh^ FWHM - relations. We want to move beyond 
single-epoch measures of the past 20 years involving random samples that yield the same 
result - quas ar black hole masses with mean Mbh"^ 8.5±1.5 (see the paradoxical result of 
Croomllionl v We think that the best path to a ccomplish this goal i nvolves generation of 
median composite spectra in the 4DE1 context (jSulentic et al.ll2002l . The composites 
can provide very high s/n line profiles from which more accurate measures can be obtained. 
Results about line profiles and trends are given in ^ Our findings are discussed in §[S]that 
presents the necessary cautions and calibration factors that can ultimately yield improved 
A/bh values. We will only brie fly analyze (ijS.SIl a lum i nosity -dependent scaling law such 
as the one recently presented in lTrakhtenbrot fc Netzeii (12012 ) as well as the determination 



of J- that is a distinct and very complex proble m (e.g., 
20 111 : iNetzer fc Marzianillioiol : IPark et all 120121 ) . 



Colhn et all 120061 : Graham et al 



2. Sample Selection and Spectral Binning 

2.1. Sample Selection 

Spectral binning offers the possibility to generate much higher s/n composite spectra us- 
ing the SDSS database if a binning context can be identified. We searched SDSS-DR7 
for sources catalogued as type-1 AGN (quasars) in the redshift range 0.4 - 0.75 and with 
magnitudes brighter than «18.5 in the g, r ox i bands. Below the assumed magnitude limit 
the low s/n of individual spectra make it difficult to estimate even FWHM H /3 and -RfbII 



well e nough to permit reliable bin assignments. We also searched the list of I Zhou et al 



(|2006l ) in order to include type-1 AGN showing broad line FWHM H/3< 1000km s'^. Such 
sources are not identified as quasars in SDSS but show strong FeH emission characteristic 
of type 1 AGN. The resultant sample consisted of 716 quasars. The next step involved dis- 
carding very noisy spectra and some sources with unusually red colors (obviously reddened 
quasars where M gllA2800 was very weak or absent) thus reducing the sample (the same 



one considered in lPaper If ) to 680 quasars. Narrow and broad absorption line sources were 
excluded from the sample unless, in the case of narrow absorption lines, the absorption 
was judged insufficient to affect our ability to measure broad line binning parameters. 

2.2. Spectral binning 

Spectral binning in the 4DE1 optical plane maps sources using measures of FWHM H/3 
vs. i?Fcii = W(FeiiA4570) / W (H/3) « I(FeiiA4570)/I(H/3), where H/3 refers to measures 
of the broad line profile of H/3 (H/3bc) corrected for narrow line H/3 emis sion as well as 



any c ontamination from [Olll]AA4959,5007; HeiiA4686, and Fell emission (jSulentic et al 



20021 ). Recent works, including the present one consider the intensity ratio rather than 
the ratio of equivalent widths. Spectral types are defined by binning in intervals of i?Fcii 
(Ai?Fcii= 0.5, from Al [i?Fcii<0.5] to A4 [1.5<i?Foii<2.0]) and FWHM (AFWHM H/3= 
4000 km s-i, from Al (FWHM H^< 4000 km s'^) and Bl to B1++ [12000 km s-^< 
FWHM H/3< 16000 km s~^]). Separation of the sample into spectral types will allow us 
to avoid the danger of mixing together sources that show very different spectral line profile 
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properties. Insofar as H/3 and Fell measures reflect physical conditions in the BLR, the 
bins isolate sources with similar values of density, ionization parameter, and L/ L-^dd (see 



Marziani et al ] |2O10l . 120011 ). Binning allows us to also take into account the complexity 



e.g. 

of line profiles and to separate, at least in an heuristic way, major flux contributions that 
are due to gas that is not likely to be virialized (i.e., the redshifted very broa d compo nent 
of H/3 in Population B sources). Bin assignments were made as described in lPaper J . We 
made a final inspection supplemented by IRAF SPLOT measures, for sources falling near 
bin boundaries. A few dozen sources were moved into adjacent bins but none of this fine 
tuning is likely to have had much effect on resultant median composites. Two bins (A2 
and Bl) are most densely populated while our binning contains enough sources to permit 
generation of composites for 8 bins that span the FWHM and i?Feii ranges occupied by 
the majority of low redshift type 1 sources ( with an inferred 2 dex r ange in Edditington 



ratio, 0.01 ^ log L/Lpdd^l that is the one of luminous type-1 AGN; iDavis fc Laoi 



2011 



Kollmeier et al 



l200d : [Marziani et aLll2003bl : IWoo fc Urrvlbooj ). 

Table [T] gives the number of sources (A'tot) in each bin. Boldface numbers indicate 
spectral types for which a composite spectrum was constructed. The separation into Pop. 
A and B source s divid es almost evenly low redshift z <1.0 quasar samples (45% Pop. B. 

(j201Gl ) and 54% in this st udy). Bin B2 was almost devoid of sou rces in 



Zamfir et al 



previous studies of low-luminosity sources ([Marziani et al 



2003a 



Zamfir et al.ll201Gl ) while 



it accounts for 3% of the present sample. Sources with B2 spectra are more numerous at 
high luminosity wher e they account for ^20% of our Very Large Telescope ISAAC sample 
( Marziani et al.l[2009[ ). Given the shape of the source distribution in the 4DE1 optical plane 
and the crudeness of the binning it is likely that bin Al is a blend of A2 and Bl sources. 
The B2 composite is mainly due to spectra of sources whose location in the 4DE1 optical 
plane is close to the boundar y of A2 and Bl. In oth er words it could be that B2 sources are 
a physically diverse mixture (jMarziani et al.ll200ir ). More than anything else, B2 provides 
a hints of the confusing results that one would obtain mixing Bl and A2 spectra. Data 
and measures on B2 are reported but scantily discussed. 



2.3. Radio-Loud Sub-Samples 



The entire sample of 680 binned quasars was cross-correlate d against the Faint Images 
of the Radio Sky at Twenty-Centimeters (FIRST) catalogue ([Becker et al.lll995r ). Radio 
maps and catalogue data for 163 FIRST-detected sources were used to isolate radio-loud 
(RL) quasars fro m the general radio-quiet (RQ) popula tion. The majority of RL sources 
belong to Pop. B ( Sulentic et al.|[2003 : Zamfir et ahlboosh . Most radio-detected RQ sources 
show core-dominated (CD), sometimes core-jet morphology. Sources were classified as RL 
lobe-dominated (LD) when a visual inspection revealed a complex, usually double-lobe 
structure. All sources except on^have radio power log Pi, > 31.6 at 1.4 GHz (P^ in units 
of erg s~^ Hz~^), the nominal lower specific power for Fanaroff-Riley II (FRII) sources. 
The radio morphology defines the classical RL phenomenon with the weakest FRII sources 
providing an empirical lower limit for this activity. Compact sources were considered core- 
domi nated (CD) RL if t heir rest frame specific emission was above the FRII lower limit. 



Zamfir et al. 



{ 2008 ). This corresponds roughlv to Rk =70 as defined in Kellermann 



et al. ([1983 ) . The lower limit of radio power in FRII sources excludes all or most radio- 
detected RQ sources that appear as CD in FIRST. Table [T] reports the number of sources 
in each bin along with the numbers of FIRST detections (Nfirst)- FIRST non-detections 
(i.e. undetected RQ sources) and detections are tabulated in columns 3 and 4. FIRST 
should detect all RL sources in the redshift range under consideration. Radio detections 



^ SDSS J093704.05+293704.9 has logP^ « 31.2 at 1.4 GHz but classical FRII morphology. 
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were assigned an FRII or CD morphology (numbers of sources are in columns 5 and 6) 
based on inspection of FIRST images. 



3. Generation and Analysis of Composite Spectra 

Three steps preceded generation of composite spectra: 1) assignment of sources to spectral 
bins as described in the previous section, 2) deredshifting of source spectra using SDSS 
z values, and 3) applicati on of a second-order c orrection connected with a known bias in 
SDSS redshift estimates (jHewett fc WildI 120101 ) . Accurate rest frame definition can be a 
serious problem. Knowledge of the rest frame is important for identification and physical 
interpretation of internal broad and narrow line shifts that show a wide velocity range . 
They differ widely from source to source and depend on ionization state ( Gaskell 19821 : 
Tvtler fc FanI Il992l ) . The technique used for redshift computation in SDSS is reliable 
but n ot accurate enough to yield uncertainties as low as « 30 km s^^ (jHewett fc Wild 
2010l ). The [Olll] AA4959,5007 line strength is a known El correlate with [Olll] AA4959,5007 
becoming weaker, broader, and hence more difficult to measure as one proceeds from B to 



the A l bin and along the sequence from bin Al to A4 (jSulentic et al 



2OO2I : IZamfir et al 



2010|). Bin A3 and A4 sources show a systematic [Olil] AA4959,5007 blueshift with respect 
to low ionization narrow lines. The shift amplitude in these "blue outliers " is also believed 



to be related to the El sequence ( Hu et al. 2008 : Komossa et al. 20081 : Marziani et al 



2003bl : lZamanov et al.ll2002 

Wavelengths of the three most prominent narrow lines [Oii]A3727, H/3, and [Oiii]A5007 
were measured whenever detected. Systematic shifts were computed taking an average of 
the three lines in each source spectrum and clipping individual measurements in cases 
of disagreement because of poor data or intrinsic blueshift of [Oiii]A5007. This means 
that in several cases only two lines could be used for redshift determination. In a few 
cases (especially in A3 and A4 bins, where oxygen lines can be very weak) only H/3nc 
could be used or, failing that, the adopted SDSS z value. The "bias" in the SDSS redshift 
determination is larger for extreme A bins but relatively small for Pop. B and bin Al 
sources. This is not surprising, since bin Bl shows stronger and sharper narrow lines that 
permit more accurate redshift determinations with the SDSS algorithm. 

After sample selection, bin assignments, and redshift correction we generated median 
composite spectra for H/3 and MgllA2800. RL and RQ composites was also constructed 
for bin Bl. Wavelength uncertainties of composite spectra were quantified by measuring 
a posteriori additional low-ionization narrow lines. Systematic effects were found to be 
small, with a radial velocity difference Aur^ 25 km s^-'^. The dispersion in rad i al vel ocity 
is CT^t,, ^ 50 km s^^ for all bins in line with estimations of lHewett fc WildI (|2010l ). In 
practice this means that any line is assumed to be unshifted if (tad,- ^ 100 km s~^. 



3. J. H/3 snd Mgu\2800 spectral range 

Interpretation of the H/3 and MgiiA2800 sp ectral ranges closely follo ws previ ous work 
(jBoroson fc Greenlll992l : iMarziani et al.ll2003ar ) and has been discussed in lPaper j . We only 
mention that emission blends near MgiiA2800 are mainly due to Fell which is modeled 
with the FeiIttt,/ emission template provided bv the last model in Table 1 in Briihweiler 



fc Verner (j2008l ) . Main results presented in this pape r were also tested with an additional 
empirical template produced bvlTsuzuki et al. Fel emission was detected in Paper 



I, but appears to be fainter ( along with any contribution due to Balmer continuum) in B 
spectral types. Adopted ranges for the multicomponent fitting analysis are 2600 A - 3050 
A and 2400-3100 A for Pop. A and B respectively. Known non-iron lines in this region 
include [Neiv] 2423.8, [Oil] 2471.0, semi-forbidden Ahi] 2669.95, O511 2672.04, and most 
probably Hei 2945.11. 
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3.1.1. The ratio of the MgiiA2800 lines 

The intensity ratio of the MgiiA2800 doublet M = /(Mg iiA2796.35) / /(MgliA2803.53j^ 
due to transitions i — Si should be 2:1 in the pure optically thin case. In the case 
of infinite optical depth the ratio should be 1:1 (i.e., if lines are fully thermalized) . For 
intermediate conditions the ratio depends on electron density, temperature, and column 
density. Under standard assumptions the condition for thermalization on the product 



ionization parameter t imes hydrogen density nu is Unn }^ 1.7 10^ cm ( Hamann et al 



19951 : iLaor et al.lll997l ). We therefore expect an R value close to 1:1 for likely values of 



U and rtH in both the BL R and the very broa d line region (VBLR) associated with the 
VBC. An array of CLOUDY ( Ferland et al. 1998h simulations predicts the intensity ratio as 
a function of hydrogen density nn and ionization parameter U (Fig. (TJ a fixed hydrogen 
column density of 10^'*cm~^ and solar abundances are assumed). 

The individual lines are too broad and hence too blended to permit reliable direct 
measurement. A-priori values must be assumed for BC and VBC (the latter only for Pop. 
B sources) and the values could be d ifferent in lieu of the different physical conditions 
expected for the BLR and the VBLR (jlVIarziani et aLlboioh . In the BLR, \ogU = -2.75, 
lognH= 12.50 yield M w 1.3:1. The observed ratio R in a typical Population A NLSyl 
source (I Zw 1) is 1.2. We assume a standard value of 1.25, noting that a small change of 
±0.05 in the doublet ratio does not affect line decomposition and FWHM measurement. 

The intensity of MgiiA2800 is difhcult to compute by photoionization codes because the 
line is primarily emitted in the partially ionized zone, where a local approximation (i.e., via 
a local mean escape probability) treatment of radiation transfer is probably inadequate. 
It is possible that the MgiiA2800 lines become fully thermalized at the extreme optical 
depth of the LIL emitting part of the BLR, justifying the assumptions of a 1:1 ratio. In the 
following we have therefore considered both cases. Mock MgiiA2800 profiles were computed 
for 1:25 and 1:1 cases assuming a wide range of line FWHM. The peak wavelengths used 
as a reference whenever the doublet is treated as unresolved, are 2800.1 A and 2799.4 A 
for the 1:1,1.25:1 cases respectively. These values are found to be almost independent of 
FWHM. We assumed that the doublet wavelength cannot be shorter than 2799.1 A, the 
value obtained for the case M =1.5:1. Values of M <; 1.5 are unlikely for the physical 
conditions expected within the BLR (Fig. [Ij. 

In the VBLR, R w 1.5 is appropriate for —1 ^ \ogU ^ —0.5 and for high density 
values (nn^ lO-'^-'^cm"'^), and only values M 1.7 are unlikely. Considering however that 
the doublet component separation is AA « 8 A <C FWHM(VBC), the VBC can be treated 
as a single line. Under this assumption the doublet ratio is indeterminate. 



4. Results 

4.1. Broad Line Profile Shapes 

Composite spectra for Pop. B sources are shown in Fig. [21 a nd Fig. [31 after continuum 



^ig. [H a n 
Paper II ). 



subtraction (Pop. A median spectra were shown in [Paper J ) . Table [51 reports FWHM, 
asymmetry index A.L, kurtosis k and centroid at fractional intensity measures for the H^S 
and MgllA2800 broad profiles (i.e. after narrow components and Fell subtraction but before 
broad/very broad component decompositi on) of both Pop. A and B spectral types. A.L, k 



and centroids are defined by Eqs. 1-3 in [Zamfir et al.[ (|2010f ). Note that A.L is defined at 
i maximum intensity and with respect to line peak while centroids are relative to source 
rest frame. These parameter values provide an empirical description of the profile without 
any model assumption. Reported uncertainties are at 2cr confidence level and have been 



In this paper vacuum wavelengths are used. Line identification includes the optical wavelength 
of traditional use. For example, we refer to [Olll]A5007 even if the reference wavelength is 5008.2 
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computed measuring the effect on the profile parameters of changes in continuum placement 
(always less than 2%, due to the extremely high s/n). Note that full profile measures for 
MgiiA2800 doublet assumed an unresolved single line. Width meas ures for MgllA2800 can 
be co nverted to single component width by subtracting 300 km s~^ ([Trakhtenbrot fc Netzer 



20121 ). Zero point wavelengths set by [Oii]A3727 and H/3 narrow component agree within 



an rms of 10 km s~^. Uncertainties were estimated by propagating the zero point rms, 
the scatter associated with SDSS wavelength calibration and the uncertainty of continuum 
placement. 

The most important results from Table [5] are the difference in line profile properties for 
Pop. A and B. FWHM differences among Pop. B spectral types arise from the definition 
of 4DE1 and how it is binned. Pop. A sources show broader, more symmetric and less 
shifted H/3 profiles than MgllA2800 with the exception of bins A3 and A4 where a blue 
asymmetry is observed. MgllA2800 and H/3 both show a symmetric unshifted profile in 
bins Al and A2. The broad MgllA2800 profile appears symmetric but blueshifted by «- 
300 km s~^ in bins A3 and A4 and these are the only sources where FWHM MgiiA2800 
exceeds FWHM H/3. This has been seen be fore and interpreted as evidence for outflows 
in the highest L/LEdd A3 and A4 sources ([Paper t ). possibly because of a distinct LIL 
component connected to the blue shift observed in HILs like ClvA1549. 

Pop. B s ources are largely represented by bins Bl and B1+. Figure [2] confirms pre- 
vious claims (jWang et al.lboogh that FWHM H/3 is systematically broader than FWHM 



MgiiA2800 by approximately 20%. This appears to be true for all B bins as well as bins 
Al and A2. Both H/3 and MgliA2800 Pop. B profiles show median red asymmetries; 
MgllA2800 show perhaps slighly weaker asymmetries. Centroid shifts are redward and of 
larger amplitude in H/3 than MgiiA2800. The red asymmetry in the lower ha lf of the profile 



was previou sly interpreted, in the case of H/3, as a signatiue of the VBC (jSulentic et al 



2002, l2000d). Apparently MgllA2800 shows a (weaker) VBC component. Much of the 



excess HP profile width may be due to a stronger VBC component. In other words, the 
measures on line profiles of Pop. B reported in Tab. [5] suggest that H/3 is more strongly 
affected than MgllA2800 by non-virial motions especially toward the line base. However, 
the VBC cannot be a full explanation for the difference in line width between MgiiA2800 
and H/3 because Al and A2 (where no VBC is present) show the same effect. 

Somewhat higher kurtosis values are found for broader Po p. B sources since Pop. A 
profiles are more sharply peaked fc.f. iKoUatschnv fc Zetz see ti5.3l for interpretation). 

The kurtosis index is consistent with a single Gaussian profile (k « 0.44) only for spectral 
type B1++ for both H/3 and MgilA2800. B1++ profiles are rare, very broad and complex 
sometimes showing double peaked lines although not in our sample. 



4.2. Multicomponent Analysis 



IRAF SPECFIT (|Krisslll994[ ) was used to model and decompose line b lends an d components 
in Pop. B composites as previously reported for Pop. A composites ( Paper I ) . We assu med 
a weak MgllA2800 narrow component in Pop. B composites (see Wang et al. 20091 and 



references therein) constrained to not exceed the ratio MgiiA2800/H/3Nc~ 2 (appropriate 
for logJ7 ~ —2 and log7iH~ 4 — 5). A double Gaussian model involving an unshifted 
BC and redshifted V BC was used for fits to both lines following previous work on H/3 
([Marziani et al.ll2009l : Fig. [2l and Fig. [3]). In contrast to Pop. B, the majority of Pop. 
A profiles were well fit by an (almost) unshifted symmetric Lorentz function. Table [3] 
presents BC and VBC measures (intensity, peak shift and FWHM) for H/3 and MgiiA2800 
for B spectral type and for the source PC 1201+436. Uncertainties are 2a confidence 
level. Uncertainties of peak shifts were computed by quadratically propagating errors on 
zero point, wavelength calibration rms and uncertainty provided by the fitting routine 
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(jPaper J ). Uncertainties reported for the VBC should be considered as formal errors and 
are computed following the assumption that that the best fit is the correct one and that 
the functional form describing the line component is also correct. The VBC is a broad 
component with poorly constrained centroid and width. More realistic uncertainties are 
estimated around 20% for both shift and width measures. 

Our operating assumption is that unobscured gas in virial motion should give rise to 
a reasonably symmetric/unshifted line profile. The EC component is therefore the virial 
estimator that we are trying to isolate. The use of a double Gaussian model in SPECFIT 
was motivated by the inflected appearance of the H/3 profile (this is evident in the Bl 
panel of Fig. [5]). Decomposition of the BC/VBC blend was the most uncertain part of the 
modeling. We derive the ratio of BC / VBC intensity for the two lines from the fit in the 
Bl, B2 and B1+ cases. However, we do not know the shapes of the individual components. 
A highly redshifted VBC is unlikely to show a symmetric Gaussian profile. Past work 
has focussed on H/3 and has used very high s/ n composites like the ones presented here . 
Past work involving BC+VBC decomposition (jMarziani et al.ll2009l : ISulentic et al.l 120021 ) 
is hmited and shows FWHM H/? BC and shift values in the range 4000-4800 km s~^ and 
zero km s^^, respectively, with the relative strength of the VBC possibly increasing with 
source luminosity. 

Modelling is clearly more ambiguous in sources with broader multicomponent profiles. 
The MgllA2800 profile in Pop. B sources also deviates from a single Gaussian show- 
ing a narrower core and more prominent wings. In other words it is more similar to a 
Lorentz/ Voigt profile that can be approximated by the sum of the Gaussians. An attempt 
to fit MgiiA2800 in the Bl composite with a single Lorentzian will reproduce the observed 
red wing but model the blue one with a significant excess. The sum of two Gaussians, with 
the broadest one slightly redshifted provides a satisfactory fit (Fig. [2]). The introduction 
of a small shift in the broader Gaussian accounts for the excess flux in the red wing, promi- 
nent in H/3 but also detected in MgliA2800. Nonetheless the assumption that the VBC is 
a shifted symmetric Gaussian is unlikely to be correct for such a redshifted component: 
non-virial motions are likely to strongly affect the VBC. We currently have no empirical 
basis for fitting the VBC with more complex models. Our best current fits suggest that 
the MgllA2800/H/3 intensity ratio differs for the VBC and BC: with MgllA2800/H/3 > 1.5 
for the BC and MgllA2800/H/3 < 1 for the VBC. This result is consistent with the idea 
that th e VBC arises in a higher-ionization VBLR near the inner edge of the BLR f Korista 



& Goadl200J). 



Modeling the VBC as a distinct component is further justified by the discovery of 
sources where the H/3 profile is dom i nated by the VBC wi th little or no detecte d BC 
(e.g. PG 1416-129 in ISulentic et al.l |2000cl ) and 3C110 in iMarziani et all |20ld )). In 



the course of assembling the sample under consideration in this study we found a few 
more VBC dominated sources. The best example involves PG 1201+436 whose H/3 and 
MgllA2800 profiles are shown in Fig. Most of the broad redshifted H/3 profile is almost 
certainly VBC with little or no BC emission. The most straightforward interpretation of 
the MgiiA2800 profile is, conversely, BC dominated with much weaker VBC emission. We 
can use these considerations to cautiously guide SPECFIT in modelling of MgiiA2800 in 
bin B sources with the broadest profiles i.e., B1++. Bl"*""*" sources however constitute less 
than 3% of our sample. They are however relevant if for no other reas on that bin Bl"'"^ 
shows the largest RL fraction in the 4DE1 sequence (jZamfir et al.ll2008l and Tab. [ij. This 
i s also the domain of some of the best known sources that show double-peaked profiles 
([Eracleous fc HalpernI 120031 ) although most do not. A BC+VBC profile with peaks near 
zero + redshift is not a double peaked profile. The fit for B1++ shown in Figure 3 assumes 
a BC-dominated MgiiA2800 profile and attempts to maximize MgliA2800 BC intensity and 
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width. It is interesting that the resuhing MgiiA2800 profile: 1) remains less broad than 
H/3bc and 2) requires a VBC component to produce an acceptable fit. 

Results detailed in Table[3]and Figures [5] +[3] indicate that in Pop. B sources MgiiA2800 
mimics H/3 by showing a VBC component. They also show that the VBC has a stronger 
effect on FWHM H/3 when measuring the full profile. In this sense FWHM MgiiA2800 
BC appears as a safer virial estimator than FWHM H/3. BC measures are preferred as 
virial estimators for both lines with a consistency correction factor taking into account the 
fact that FWHM MgiiA2800 (BC) is 20% narrower than FWHM H/3 for the majority of 
sources (B bins and bins A1+A2, 90 % of our sample). Nonetheless, all measures reported 
in Tab. [Hand Table |3] for all spectral types excluding A3 and A4 show that the FWHM of 
any MgiiA2800 component is narrower than the equivalent H/3 component. This result is 
robust since it also holds for the full profile. 

Tab. [3] indicates that Bl sources may show a significant MgllA2800 BC blueshift, with 
smaller amplitude in km s~^ than the one measured in Pop. A. A small blueshift might 
also be present for B1++. Part of the BC blueshift might result from modeling the VBC 
with a symmetric profile and possibly over-subtracting flux from the red side of the BC. It 
is however unlikely that this accounts for the full effect since a blueshift is also measured 
for the full profile: a blue-side excess is visible in the median Bl MgllA2800 profile shown 
in Fig. m Tab. [2] indicates a MgiiA2800 blueshift for spectral types Bl and B2 with 
any blueshift in Bl''"*' BC not detected by our parameterization of the full profile. The 
relevance of the blue shifts in Pop. B profiles will be further discussed in t)4.4l and ^ 

4.3. Radio Loud 

RL sources, both CD and FRII (Fig. [5]) Bl medians show a symmetric (FRII) or redshifted 
(CD) MgiiA2800 line core. The profiles resembles the typical Bl RQ and full sample profiles 
but there is no hint of possible blue shifts. A shift to the red is observed at 1/4 intensity 
level but, again, the shift amplitude is less than for H/3. The redward asymmetry is 
especially prominent in CD sources with a ^ centroid displacement reaching 600 km s~^. 
CD and FRII sources do not show any large difference in FWHM since they were both 
chosen to belong to spectral type Bl by imposing the same condition on FWHM H/3. Their 
interpretation will be discussed in ti5.4l 

4.4. Major trends 

We focus on line shifts and widths because they are parameters relevant to the credi- 
bility/definition of a virial broadening estimator. Comparison of centroid shift values at 
different heights in a line give the best way to evaluate asymmetry and line shift. The 
behavior of Pop. B is different than Pop. A showing a redshift that is strongest in the 
line base resulting in red asymmetric profiles. We assume that modeling and subtraction 
of a VBC component will enable us to measure the width of a symmetric and unshifted 
BC virial estimator. Pop A sources in bins Al and A2 yield a FWHM measure requiring 
no VBC correction. A correction is needed only to make consistent virial estimate us- 
ing FWHM H/3 and FWHM MgiiA2800. In bins A3 and A4 the FWHM H/3 measure is 
preferred after correction for a blue asymmetry. 

Rather than the amplitude in km s~^, a value normalized by line width provides a more 
direct evaluation of the "dynamical relevance" of any line shift. We can define the following 
parameters 





,i = 1,2,3, 



(2) 



FW|M 
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where the centroid c(^) at fractional peak intensity j has been normahzed by the full 
width at the same fractional peak intensity FW|M. Fig. Oshows the behavior of (5(|) and 
of 6{j) as a function of spectral type (error bars are obtained propagating quadratically 
the uncertainties in Table. [5]). A radial velocity v^- shift close to a line peak may represent 
a systematic effect involving the entire emitting region only if it is also seen at ^ and 
J fractional intensities. Fig. [H] confirms that the Bl MgiiA2800 blueshift is much less 
relevant than shifts found in bins A3 and A4. We also note that in bin A3 and A4 the 
shift measured for the full MgliA2800 profile is probably a lower limit, since the profile can 
be interpreted as the sum of an unshifted line + a blueshifted component (^3.3 of Paper 
I). The line base might be affected by the presence of additional kinematic components 
that, in our interpre tation of the line profiles, are the blueshifted component ( blue in 



Marziani et al 



and the VBC for Pop. A and B respectively. Considering the j 



centroids as a function of spectral type, only bins A3 and A4 show a relatively large H/3 
blueshift (the peaks are almost unshifted). Bin Al and A2 show symmetric profiles in 
both H/3 and MgilA2800. The situation changes in Pop. B where the Bl and B1+ H/3 
profiles are strongly red asymmetric with the B1+ centroid shift reaching « 800 km s~^. 
Conversely the MgiiA2800 profile shows a more symmetric shape. There is a significant 
shift at J intensity level but its amplitude is ~0.25 that of H/3. Given the large FW-jM for 
MgliA2800, the 5 values are lower than in the case of H/3. Similar considerations apply to 
i intensity - the line intensity level at which the virial broadening is computed: MgiiA2800 
shows modest amplitude normalized shifts (Fig. [7]). 

Eddington ratio estimates are needed to infer the physical meaning of differences be- 
tween A and B spectral types. We computed the median mass from the medians of the 
fiuxes of individual sources in each bin and from FWH M m easured on the med ian compos- 



ite. We follow the prescriptions of I Asset et al.l (|201lf ) and lShen fc Liul (|2012l ) to compute 



Mbh from H/3 and MgiiA2800 line widths respectively. Trends are preserved also with older 
relations linking the black hole mass to the FWHM and the continuum at 3100 or 5100 
A. The bolometri c luminosity needed for L/LEdd has been derived from the luminosity at 
5100 A following ETemmen fc BrothertonI |201oh . Values of bolometric luminosity, A/bh 



and L/LEdd computed using the BC as a virial broadening estimator are reported in the 
first columns Tabled followed by Mbh and L/L^dd computations using the whole profile 
(treating MgiiA2800 as a single line). Mbh values derived from BC and presented in Table 
m show a trend toward decreasing A/bh from Al to A4. A constant median L in the bins 
(logi ^ 46.0) implies a trend of increasing L/L^dd from Al to A4. We do not see this 
trend if we use MgiiA2800 measures reflecting the absence of a clear correlation between 
H/3 and MgliA2800 line profile measures in Pop. A. Pop. B sources show a more regular 
behavior preserving a trend. This is mainly due to the definition of B spectral types based 
on increasing line width. However, inclusion of VBC leads to an overestimate of log Mbh 
for both H/3 and MgllA2800 derived masses. 

Fig. [7]shows the behavior of S{^) and of the FWHM ratio between the broad component 
of MgilA2800 and H^ as a function of Eddington ratio since the El sequence is mainly a 
sequence of £/LEddlj The parameter S{^) is considered as an appropriate indicator of the 
dynamical relevance of the shift because FWHM is a virial broadening estimator of choice. 
The S trends at j and | peak intensity as a function of spectral type are preserved also as 
a function of L/ LEdd and are not shown again. The MgllA2800 blueshift starts increasing 
at log L/L-Edd^ —0.6, as the A3 and A4 sources behave differently from all other sources. 



Our sample bins have median Lboi that is almost constant. This follows from the flux and 
redshift limits of our sample, and implies that for any trend as function of L/LEdd, there is a 
corresponding trend with Mbh- A/bh is a one-way estimator of a quasar evolutionary status. 
However, the presence of line shifts is more probably associated to the relative balance of grav- 
itation and radiation forces so that we will assume in the following that the relevant physical 
parameter is L/LEdd- 
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Paper ll . The marginal blueshift detected at | maximum for spectral type 



as discussed in 

Bl becomes untraceable at ^ intensity. At ^ maximum the MgiiA2800 profile is mostly 
unperturbed in all B bins (save B2 that includes a tiny minority of sources and is probably 
a confusing mixup), while the H/? profile shows a highly significant shift to the red. 



5. Discussion 



5.1. Mgu\2800 and H/3 in the Pop. A and B Context 

Virial estimat es of Mbh for large quasa r samples have been avai lable now for more than 
10 years (e.g. iMcLure fc Dunlod l2004l : iMcLure fc Jarvid 120021 ) . Estimated masses lie 
within the range log Mbh^ 6.0 - 11.0 with uncertainties i.n the range 0.2-0.4 dex at la 
confidence. The latter uncertainty value appears more realistic when different estimates 
for the same sources are compared. This is also confirmed when we compare different mean 
mass estimates for quasar populations (e.g. radio-quiet vs. radio-loud). This situation has 
led to claims that FWHM measures contain little or no information about black hole mass 



(|Croomll201ir ) an d others arguing that FW HM measures for different lines are equally valid 
Mbh estimators (jVestergaard et al.ll2011l ). We argue that neither of these claims is valid 
and that the path to a clearer picture lies within the concept of a quasar parameter space 
like 4DE1. Quasars do not distribute randomly in 4DE1. 

How can we move towards more accurate and consistent Mbh estimates for individual 
quasars and quasar populations? It is unl ikely that incons istencies are due only to com- 
parisons involving spectra of different s/n (jAssef et al.ll2012l ). Inconsistencies remain when 
comparing estimates using the same (i.e. SDSS) spectra. The s/n is certainly an issue but 
its effect depends upo n the structure of the line profiles used as virial estimators (see e.g. 
Shen et al. 20081 . 2011 ) . Lines well-fit with single unshifted components are likely the most 
robust in the face of declining s/n. However many broad lines in type 1 AGN show line 
shifts and asymmetries while others are obviously not well fit with a single component. 
Fitting lines with a single Gaussian is inappropriate because the profile shapes are rarely 
well approximated by a single Gaussian. For example, an attempt to fit the MgliA2800 
line of spectral type Bl with a single Gaussian would lead to a (bad) fit resulting in an 
overestimation of FWHM by 700 km s^^. Effects of erroneous assumptions on fitting func- 
tion depend on s/n in a way that depends on the intrinsic shape and on noise properties. 
Composites spectra have very high s/n and allow one to focus on systematic effects with 
minimal measurement (statistical) errors. The challenge is to generate composites in a 
clear and unbiased way which means binning physically similar quasars. 

If any broad lines serve as virial estimators then MgllA2800 and H/3 are the safest 
choices. Results presented in this paper suggest that MgiiA2800 is the line least affected 
by profile shifts and asymmetries in 90% of all quasars. We maintain that a first step in the 
path to progress requires recognition that quasars showing FWHM H/3 < or > 4000km s~^ 
(Population A and B respectively) have very different profile properties! Collin et al. 20061 : 
Sulentic et allboOObl : IZamfir et allboioh . In other words, the dispersion in estimated Mbh 
and resultant LEdd measures is not found by luminosity binning (at fixed z) but rather by 
binning in the 4DE1 context that reflects differences in BLR kinematics and geometry. 

Figure [5] (left) overlays our composite measures of FWHM H/3 vs. FWHM MgiiA2800 
(with the MgiiA2800 dou blet treated as a single line) on the individual source measures 



from I Wang et al.l (12008 ). Full profile measures follow the source distribution and the 
best- fit regression of IWang et al.l (|2009f ). The large scatter at the low FWHM end and 
the convergence of the regression towards parity refiect the unusual behavior of sources in 
bins A3 and A4 where the FWHM offset between H/3 and MgiiA2800 breaks down. The 
convergence is also affected at the high FWHM end where the prominence of the VBC 
component increases in H/3 more than in MgllA2800. Figure \E\ (right panel) presents our 
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BC values (all Pop. B sources corrected for the VBC) along with the parity line, the Wang 



et al. regression and our best fit to the median composite values. The main result 

is that we find no convergence but rather a constant offset (MgiiA2800 «20% narrower) 
from bin A2 through all B bins. 



5.2. Implications for Mbh estimates of quasars: a tentative recipe 
The previous analysis suggests the following recipe. 

1 Identify a source as Population A or B - on the basis of FWHM H/3 < or > 4000 km s~^ 
respectively. 

2a Bin A3 and A4 MgilA2800 profiles are unsuitable as virial estimators. Bin A3 and 
A4 H/3 profiles are suitable if corrected for strong Fell contamination and for a blue 
component. 

2b Al and A2. Both and MgllA2800 can be used as estimators without any correction. 
3 Pop. B. Apply correction because of VBC to both H/3 and MgilA2800, as described 
below. 

A correction factor can be defined as the ratio between observed FWHM and BC 
FWHM: 



FWHM(line),i,. = eFWHM(hne),i,, (3) 

where 
^ FWHM(hne)Bc 

FWHM(line)„^,' ^ ' 

and FWHM(line)Q|^j, is measured on the whole profile (treating the MgllA2800 doublet as 
a single line). If a single-epoch FWHM measurement is made without consideration of 
the intrinsic shape then ^ is as reported in Tab. [S] values are derived from Tab. [3] and 
Tab. [5]). The ^ factor is less important for MgllA2800 because the VBC component is 
always weaker. In the case of MgllA2800 a suitable ^ value is « 0.85. Table [5] indicates 
that a simple correction factor ^ « 0.75 — 0.8 could be applied to H/3 once the object is 
classified as Pop. B, as suggested by Fig. |8l A simulation degrading the Bl spectrum with 
a combination of Gaussian and Poissonian noise shows that it is p ossible to recover a basi c 



estimate of FWHM (with heavy smoothing) down to s/n « 3 fc.f. IShen et al.ll2008l . l201ll ). 
even if information on profile shape is completely lost. We also note that the f correction 
may be accurate excluding t he most luminous sou rces where the VBC is extremely strong 



i.e., for iboi^ lO'^^ergs s ^ ( Marziani et al. 20091 ) . or the profile is not too irregular as it 



is the case of several sources in B1++. The third column of Tab. [S] lists the multiplicative 
factor needed to convert from observed FWHM H^ to FWHM MgiiA2800 BC. The last 
column lists the factor needed to convert the observed MgiiA2800 FWHM into FWHM 
H/3bc. 

The recipe reported above accounts for our result that the difference between MgllA2800 
and H/3 Mbh is positive at large L/LEdd, and negative toward the lower end of the L/Z/Edd 
distribution. If no correction is applied then MgilA2800 will be measured broader than H/3 
toward high L/^Edd values (20% in A3 and « 50% in bin A4). Therefore at high L/LEdd 
the MgiiA2800 masses will be overestimated with respect to H/3 which should actually 
provide a more reliable estimate. At low L/L^dd the effect of the VBC on FWHM H/3 
is stronger yi elding an MeHestimate lar ger than the one from MgiiA2800. The agreement 
with Fig. 1 of Onken fc KoUmeier ( 2008h is quantitative (if rare outlying points in their Fig. 



1 are excluded; medians by definition give no or low weight to extremes in a distribution) : 
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A log AfBH= log AfBH(MgiiA2800) - log Mbh(H/3) w 0.4 dex for bin A4 with logL/iEdd- 
-0.2, and A log Mrh- -0.4 at logL/LEdd« -1-7. 

Even if A3 and A4 represent only wlO% of our low redshift sample, high L/^Edd are 
discovered with increasing frequency in high-z, flux limited samples. A recent study for 
2; <1.8 sources involving MgiiA280G and Ha finds a small offset in the A'/bh estimates, in 
the sense that MgiiA2 800-derived masses are slightly larger with respect to those from H/3 
([Matsuoka et al.ll2013r ). This might indicate an increasing frequency of MgllA2800 profiles 
affected by non-virial, outflow motions. 



5.3. Origin of MgiiX2800 emission and profile blue shifts 

A major result of this paper is that FWHM MgllA2800 is systematically narrower than 
FWHM H/3, excluding types A3 and A4. This holds for the full profile as well as for all the 
components identified in the lines: i.e. BC, VBC, blueshifted component. The simplest 
explanation is that MgiiA2800 is emitted at larger radius, on average, than H/3 assuming 
that the gas is photoionized by a central continuum source. The emissivity S of the two 
lines was computed as a function of the ionization parameter U using cloudy simulations 
described in ti3.1.1l The dependence on U was then converted to radial distance assuming 
a constant column density log riH= 11 [cm^'^]. If we approximate S with a power-law then 
S oc yields i;(MgllA2800) oc r-°-'^ with a steeper trend for H/3, S(H/3) oc r^O-^, if 
logU ^ —0.5. Elementary considerations suggest that the different emissivity law can lead 
to a different shape in the line wings and hence to a different FWHM for the two lines. 
The energy emitted in the line at distance r can be written as E{r)dr — 4:TTr^J^{r)fc{r)dr, 
where /c is the covering factor. Since the virial assumption implies that r (x the 



specific energy emitted per unit radial velocity is E{v)dv oc v' 



^dv 



if /c(r) oc r"^"*"'. For a = 0.5 and q = 0, E{v) oc which is the same general power- 
law describing Lorentzian line wings. A flatter emissivity law will yield a more peaked 
profile. Using the emissivity laws for Mgll A2800 and H /3, a 20% width ratio is obtained for 

199nh . A 



c.f. Netzer 



reasonable values of(7(— l^gSji, , 

the computation of the line profile following a weak- field approximation ( Chen et al. 19891 : 
Sulentic et al. 19981 ). Using the CLOUDY-computed emissivity laws for H/3 and MgllA2800, 
the profile of MgllA2800 turns out to be systematically narrower than H/3 by ^10%. The 
width ratio can be «20% if the inner radius of the emitting disk is k, 250 gravitational 
radii. Therefore, the line width ratio and the Lorentzian-like profile shape can be explained 
in the context of virialialized gas motions. 

The MgiiA2800 blueshift is most straightforwardly attributed to outfiow of the l ine 
emitting gas with preferential obscuration of the receding part of the flow iPaper 1 A 
blueshift is suspected in the Bl composite (Fig. [2]) but is not significant at a 2a confidence 
level. The small amplitude of the shifts does not rule out a change in the intrinsic multiplet 
ratio as an alternative explanation in Pop. B sources: a value R = 2 would yield an 

effective wavelength close to the one of the ^Ps — Si transition, 2796.35. If the effective 

22 

wavelength is close to this value, then any systematic blueshift ;S 300km s^^ will become 
insignificant. If the blueshift is associated with a low-density high-ionization component 
then a value as large as R wl.7 could be possible. However, the absence of any blueshift 
in bin Al that is the simplest profile to analyze (1 component, weak Fell) argues against 
this suggestion. 

Outflows driven by line or ionizing photon pressure can accelerate the line-emitting 
gas to a terminal velocity that is inversely proportional to the square root of the column 
density Nc and directly proportional to the Eddington ratio i.e., 



vt oc 



Edd 



t^Kopl 



(5) 
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where UKepi is the Keplerian velocity at the outflow starting radius. 

If real the modest shift observed in the median spectra of Bl and B1++ could be 
interpreted as the terminal velocity of a radiation driven outflow. This shift could occur 
even if L/L-^dd^ 0.01 in gas with moderate column density N^^ 10^'^cm^^ (provided that 
the gas remains optically thick to the ionising continuum) that may in turn support the idea 
that resonant lines are driving, at least in part, the acceleration of the outflow ( Gangulv 
et al. 



200 



I 



Proga et al 



2001 



I 



Sulentic et al 



20071 ). There is an intriguing similarity in 



the MgllA2800 and ClvA1549 behavior that may be related to the resonant nature of both 
lines but the effects on MgllA2800 appear by far less relevant than on CivA1549 since the 
median MgliA2800 profile retains an overall symmetric appearance. The Bl and B1++ 
profiles are on average only slightly perturbed by blueshifted emission: low-A^c gas might 
produce some scatter among B sources since it is most strongly affecte d by radiative forces 
and m ay become unbound following a continuum luminosity increase (jNetzer fc Marziani 
201GI ). B1++ sources are likely to be some of the lowest L/ L-^dd sources. They are unstable 



in conti nuum emission and often show strongly variable emission line profiles with irregular 
shapes ( Lewis et al. 201ol ). They are also rare 3% in the present sample) and are 
interpreted as the most evolved, perhaps dying or starving quasars in the 4DE1 sequence. 
An example involves the well-studied FRII source 3C 390.3 which belongs to bin B1++ with 
line profiles characterized by a prominent blueshifted peak that produces a significant net 
profile blueshift. Another example with extreme velocity offset and large FWHM involves 
SDSS 0956+5128 Jsteinhardt et al.ll2012h. E ven if the inte rpretation of 3C 390.3 is currently 



a subject of debate (e.g 



Dietrich et al 



2012lMZhang||20111 ). an application of the method of 



Negrete et aP |2012h indicates relatively low-density and high-ionization consistent with the 



overall properties of th e blueshifted component identified in the spectra of many quasars 



(jMarziani et al.ll201Gl ). A significant fraction of B1++ sources in the present sample 
indeed show MgilA2800 profiles with large blueshifts/blueward asymmetries. Therefore, 
the marginal blueshifts of spectral type B1++ median composite may reflect the somewhat 
erratic nature of low-density outflowing gas detected in sources whose unstable accretion 
rate may induce significant continuum fluctuations. 

The high s/n of our composite spectra reveals a semi-broad component in the [Oili] AA4959,5007 
profiles and on the blue side of H/3nc- Bl offers the most striking case along with PG 
1201+436. In the Bl H/3 composite spectrum the excess flux is explained as a semi-broad 
component with width and shift similar to the [Oili] AA4959,5007 semi-broad component. 
That feature does not belong to the broad profile and would induce a blueshift of ^ —100 
km s~^ into the BC model or into the measures of the broad line profile if not adequately 
taken into account. A similar interpretation is also possible for the small blueshifted 
"bump" observed in the H/3 and MgiiA2800 profiles of PG 1201+436 (blue fine in Fig. [J). 
The bump is modeled with a Gaussian component of width 1500 ~ 2500 km s~^ and shift 

— 1500 2000 km s^^. The inferred velocities and widths appear rather extreme to other 

cases where a semi-broad component of [Olll] AA4959,5007 is isolated. We can speculate 
that also this source is showing outflow largely in a region that might be at the bound- 
ary between the BLR and the narrow- line region (jZamanov et al.l 120021 ). A more reflned 
analysis of PG 1201+436 is deferred to later work (Sulentic et al., in preparation). 



5.4. Radio quiet and radio loud 



Core-dominated and lobe-dominated sources show a fairly symmetric profile, with evidence 
of a reds hifted base due to the VBC that is especiallv prominent in CD sources (e.g.. Punslv 



& Zhang l201lh . The red excess is clearl y detected als o in MgllA2800 (a very convincing 



case is provided by the blazar 3C 279, iPunslvl l2013f ) although the A.I. is significantly 



lower than H/3. The absence of any hint of blueshifts in RL profiles could be related to the 
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full suppression of an accretion disk wind, or it could be due to orientation. Bl bin was 
defined on a limited range of line width. Therefore we cannot expect to see the systematic 
FWHM difference mentioned above. In addition the fraction of FRII increases toward bin 
B1++ i.e., broader sources. In bin Bl we are probably considering a subsample of FRII 
that show narrower profile then t he majority of the population. If the FWHM is affected 
by orientation as expected (e.g., 



2OO3I : IWills fc Brownelll986l : 



Rokaki et al 



Zamfiretal 



tne p 
200a 



Runnoe et al 



2012; 



Sulentic et al 



20081 ). we are observing sources that tend to be 



seen at smaller viewing angle than t he general FRII population. If the wind flow follows 
a bent path over the accretion disk ( Elvis! 200oh or if the flow is constrained closer to the 
disk plane, the absence of a blueshift might be therefore related more to orientation than 
to full suppression of disk wi nds in RL sources by t he pressure of a cocoon associated 
with the relativistic ejections (jNorman fc Milevlll984^ . Disk winds are probably not fully 
suppressed in radio-loud sources, as i ndicated by sources h ke 3C390.3 and by the discovery 
of fast outflows in radio- loud AGNs (jTombesi et al.ll2010r ). 



5.5. MgiiX2800 virial broadening and virial product 

The FWHM difference between 11/3 and MgllA2800 is very important. From the data pre- 
sented in this paper, with the exception of spectral types A3 and A4, the FWHM(H/3) ^ 
FWHM (MgiiA2800) both if the whole profile is considered or if only the broad compo- 
nent is isolated. The FWHM of a line by itself is of limited usefulness in computing Mbh 
following Eq. [1] A virial broadening estimator (let it be a line FWHM or velocity disper- 
sion) should be associated to a typical distance to put in use the simple virial relation that 
is customarily employed in Mbh calculations ( McLure fc Dunlopl l2004l : IShen et al.l 12008 ; 
Trakhtenbrot fc Netzer 2OI2I : Vestergaard fc Peterson! 2006 '). In principle, there are two 
main possibilities in the interpretation of the width differences (not mutually exclusive). 
The smaller MgilA2800 FWHM may imply that (1) MgiiA2800 is emitted farther out from 
the central continuum source according to the virial law Av oc r~^/^. If the FWHM ratio is 
« 0.8, then MgiiA2800 should occur at a distance significantly larger than ¥1/3, by a factor 
1.4; (2) only part of the gas emitting H/3 is emitting MgiiA2800. If (1) applies the scaling 
relations involving MgllA2800 should be defined with some care since tblr defined from re- 
verberation mapping of H/3 would suggest smaller distances. In the second case, MgllA2800 
can be interpreted as perhaps the "best" estimator, associated to gas very optically thick 
that responds to continuum changes. Reverberation mapping studies of MgllA2800 are not 
yet conclusive, and they leave o pen the possibility that the average emitting distanc e of 



MgliA2800 is larger than for H/3 |Metzroth et al.ll20od ; E'revese et al.ll2007l ; IWoolf2008h . In 



addition to the theoretical considerations of ti5.3l that favor a larger distance for MgllA2800 
emission, there is empirical evidence provided by the difference between the MgllA2800 
and H/3 profile shape. For Pop. B sources, the H/3 VBC gives rise to a stronger redward 
asymmetry than in MgiiA2800, and suggests that MgiiA2800 emission is occurring farther 
away from the central continuum source. This is consistent with the observation of smaller 
FWHM(MgiiA2800) for both BC and whole profile for spectral type Al and A2 that do not 
show appreciable VBC emission. The implication of a larger emissivity-weighted distance 
for MgiiA2800 is a shift in the correlation between tblr (comp uted from H/3) and the 
lumin osity at 3000 A, by « 0. 16 dex in tblr hi Eq. 7 and Fig . 4 of iTrakhtenbrot fc Netzei 



(|2012l ). However, Eq. 12 of ITrakhtenbrot fc Netzei] 120121 ) will be unaffected since the 



virial factor is scaled to yield the A/bh derived from H/3. The difference is that T should 
be « 1.0 and not «1.3 - 1.4 that follows from the assumption that H/3 and MgiiA2800 
are emitted at the same distance. Scalings adopted in all previous studies assumed that 
a single J- value is appropriate for all quasars. This assumption is most likely incorrect, 
especially if non-gravitational forces are at play. Previous work suggested a significant T 
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difference between Pop. A and B (jCoifin et afJl2006f ). A much needed improvement is 
therefore an evaluation of as a function of spectral type or L/L-Edd- This is however 
beyond the scope of the present paper. 



6. Conclusion 

Both H/3 and MgllA2800 appear to be suitable for 7\/bh virial estimates of black hole 
mass at least for a large fraction of quasars. Both low ionization lines, however, show 
profile shifts and asymmetries that must be taken into account if our goal is to improve 
the poor present-day accuracy of most estimates. Right now, this analysis suggests that 
masses for Pop. B sources based on ¥1/3 are likely systematically overestimated. Further 
progress therefore requires subdivision of large quasar samples into Pop. A and B with 
further subdivision into spectral types facilitating consideration of Mbh and/or L/L^dd 
trends. 4D Eigenvector 1 currently offers the most effective context for binning. Such 
source discrimination cannot be accomplished with L, z binning because sources with very 
different line profile properties show similar Lboi even when they host super-massive black 
holes with significantly different masses (i.e. the highest accretors are connected with 
Pop. A sources involving the smallest black hole masses-effectively blurring any expected 
mass-luminosity correlation for quasars). 

Our ability to distinguish between Pop. A and B sources is fortunately preserved even at 
low s/n. We defined simple prescriptions with the goal of providing a better approximation 
of the virial broadening for single-epoch observations. 
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Table 1. Sample: number of sources in each spectral type 



Sp. TypeA^tot 




/\rb 

-"^FIRST 


Arc 

^^FRII 


^CD 


Al 


97 


79 


24 


10 


8 


A2 


156 


152 


20 


1 


3 


A3 


43 


39 


14 





4 


A4 


15 


13 


4 





2 


Bl 


218 


179 


58 


23 


16 


B1+ 


115 


95 


27 


11 


9 


B1++ 


17 


9 


11 


4 


4 


B2 


19 


16 


5 





3 


All 


680 


582 


163 


49 


49 



^ Number of radio-quiet sources defined as -/Vrq= A^tot - ^cd - ^frii- 

^ Number of FIRST-detected sources. Note that iVpiRST > A^^cd + A^frii- 

Number of FIRST-detected sources with compact morphology and emitted rest frame 
power logP,y > 31.6 {Pi, in units of erg s~^ Hz~^). 

Number of FIRST-detected sources with extended morphology. 
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Table 2. Measured Parameters on the H/3 and MgllA2800 Profiles of Composite Spectra 



Sp. Type H/3 MgiiA2800 





FWHM'^ 


A.I.'' 










c(0.9)^''^ 


FWHNP 


A.I.'^ 




,(r)a,a 






c(0.9)'''^ 














Pop. 


A Full Sample 














Al 


3180±90 


0.00±0.04 


0.33± 0.02 


50±110 


50±60 


50±40 


110±50 


3040±80 


0.01±0.04 


0.36±0.02 


80±100 


80±50 


70±40 


60±50 


A2 


2900±210 


O.OOiO.lO 


0.33± 0.04 


-20±240 


-20±110 


-20±90 


-40±110 


2600±160 


0.01±0.09 


0.38±0.04 


-60±200 


-70±90 


-70±80 


-80±100 


A3 


2510±200 


-0.15 ±0.10 


0.30± 0.04 


-330±230 


-90±110 


-10±80 


30±90 


2560±160 


0.01±0.09 


0.38±0.04 


-140±200 


-140±90 


-150±80 


-160±100 


A4 


2530±250 


-0.23±0.10 


0.27± 0.04 


-530±250 


-150±130 


10±80 


80±100 


2970±210 


0.01±0.09 


0.35±0.04 


-240±220 


-250±120 


-250±90 


-260±130 














Pop. 


B Full Sample 














Bl 


7050±180 


0.09±0.04 


0.38±0.02 


660±230 


220±100 


140±80 


240±110 


5300±120 


0.07±0.03 


0.42±0.02 


230±140 


10±70 


-40±70 


-60±90 


B1+ 


9420±220 


0.08±0.03 


0.40±0.02 


790±220 


380±120 


250±110 


430±150 


7480±170 


0.06±0.03 


0.41±0.02 


380±180 


140±100 


70±90 


50±120 


B1++ 


13080±720 


0.04±0.06 


0.43±0.04 


600±570 


400±360 


290±350 


480±480 


9360±490 


0.02±0.06 


0.45±0.04 


100±360 


30±240 


0±250 


-10±350 


B2 


5470±340 


0.15±0.14 


0.35±0.06 


680±650 


60±180 


-30±160 


-100±200 


3950±230 


O.lOiO.ll 


0.39±0.05 


120±360 


-130±130 


-180±120 


-200±150 














Radio Loud 
















Bl CD 


7060±460 


0.16±0.04 


0.36±0.02 


1230±230 


440±230 


310±80 


500±240 


5690±330 


0.10±0.03 


0.39±0.02 


590±160 


270±170 


180±70 


160±210 


Bl FRII 


6790±410 


0.09±0.05 


0.37±0.02 


850±260 


420±200 


330±80 


620±250 


5500±300 


0.05±0.03 


0.43±0.02 


230±120 


100±160 


50±80 


30±210 



^ In units of km s ^ . 

Dimensionless asymmetry index computed at j peak intensity. 
^ Dimensionless kurtosis index defined as the ratio between the full width at | and ^ maximum. 

Centroids at i, i, | and maximum intensity. 



Table 3. Derived Quantities from the H/? and MgllA2800 Profiles Composite Spectra Multicomponent Analysis 



Sp. Type 












MgnA2800 




EC" 


Shift 


FWHM'^ 


VBC 


Snitt 1 WHM BC 


Shift'' FWHM^ 


VBC^ 


Shift'' 


1 WHM 












Pop. B Full Sample 










Bl 


46 


20^= 


5720±100 


51 


1740±15015600±150 63 


-165±504620±100 


36 


2160^= 


11500±300 


B1+ 


46 




7750±200 


45 


1920±19015700±260 68 


-60±80 6330±100 


51 


1520^= 


14500±1230 


B1++ 


28 


-165'' 


10100±420 


65 


915±70 15830±860 82 


-105'= 8870±420 


20 


1305^= 


13000^= 


B2 


33 


-120±60 4410±100 


40 


1950±280 13400±600 43 


-240±40 3220±240 


38 


1420±540 


10000^= 


PG 1201+436 42^ 




7500'= 


179"^ 


2125±37011800±770 118^ 


0^= 7220±590 


M"^ 


1790±330 


9780±740 












Radio Loud 










Bl CD 


34 




5330±200 


50 


1935±19013160±260 60 


45±60 4680±120 


51 


1630^= 


11000^= 


Bl FRII 


47 


235±55 


5370±130 


60 


1645±17015450±530 74 


-55±70 4790±200 


37 


1630^= 


11000^= 



^ Line intensity normalized to the continuum at 5080 A. The value roughly corresponds to the rest-frame equivalent width in A for H/3. 
within each spectral type can be computed from the reported values. 

^ In units of km s~^. 

" Kept fixed or allowed to vary within a narrow interval in SPECFIT analysis. 
In units of 10~^^ erg s~^ cm~^ A~^ 



Table 4. Median bolometric luminosity, black hole mass and Eddington ratio 



Sp. Type 






EC only 








Whole profile 






T a 
^bol 




Mbh" 












L/L^dd'^ 






K/S 


MgiiA2800 


H/3 


MgiiA2800 


H/3 


MgiiA2800 




MgiiA2800 










Pop. 


A Full Sample 










Al 


46.18 


8.67 


8.62 


-0.60 


-0.56 


8.67 


8.71 


-0.60 


-0.64 


A2 


46.16 


8.57 


8.45 


-0.53 


-0.40 


8.57 


8.57 


-0.52 


-0.52 


A3 


46.15 


8.33 


8.41 


-0.28 


-0.37 


8.45 


8.53 


-0.41 


-0.49 


A4 


46.20 


8.27 


8.56 


-0.18 


-0.47 


8.48 


8.66 


-0.39 


-0.57 










Pop. 


B Full Sample 










Bl 


46.13 


9.15 


8.98 


-1.13 


-0.96 


9.33 


9.08 


-1.31 


-1.06 


B1+ 


45.96 


9.29 


9.06 


-1.44 


-1.21 


9.46 


9.18 


-1.61 


-1.33 


B1++ 


45.96 


9.53 


9.14 


-1.68 


-1.29 


9.75 


9.35 


-1.90 


-1.50 


B2 


46.23 


8.99 


8.74 


-0.87 


-0.62 


9.17 


8.89 


-1.05 


-0.77 



^ Decimal logarithm of bolo metric luminosity in un i ts of e rgs s ^. A bolometric cor- 
rection computed following iNemmen fc Brothertonl (|2010[ ) was applied to the 5100 A 
specific luminosity. 

Decimal l ogarithm of black h ole mass in solar units. Mbh has been estimated from the 
relation of IShen fc Liul (|2012[ ) for MgiiA2800 and from the relation lAssef et all (|2011h 
for H/3. 

^ Decimal logarithm of Eddington ratio. 
Table 5. Conversion factors 







e(MgiiA2800)^ 


MgiiA2800Bc 


/ MgiiA2800 obs/ 








H/3 obs 


H/3bc 


A2 


1.00 


1.00 


0.76 


0.90 


Al 


1.00 


1.00 


0.85 


0.96 


Bl 


0.81 


0.87 


0.66 


0.93 


Bl CD 


0.75 


0.82 


0.66 


1.07 


Bl FRII 


0.79 


0.87 


0.71 


1.02 


B1+ 


0.82 


0.85 


0.67 


0.97 


B1++ 


0.77 


0.95 


0.68 


0.93 



Conversion factor ^ is the FWHM of BC divided by the FWHM of the full ("observed") 
profile. The MgiiA2800 doublet is treated as a single feature. 
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Fig. 2. Spectra of H/3 (left panels) and of MgllA2800 (right panels) for spectral types Bl 
to Bl^+ (from top to bottom). The horizontal scales are in rest frame wavelength [A] or radial 
velocity, with the origin set at the rest frame wavelength. The black lines shows the original, 
continuum-subtracted spectrum, while the dashed magenta line the model with all emission line 
components. The thick black line is the EC; the thick red line shows the VBC. The green lines 
trace the Fe llopt and Feiiuv contribution, and the gold lines various contributions associated 
to the NLR (H/3nc, [Oiii] AA4959,5007, and MgllA2800 narrow component when appropriate). 
Excess emission at A 2900A visible in this and in the next two Figures is likely due 
to a combination of Fel and Balmer continuum emission. 
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Fig. 3. Spectra of H/3 (left) and of MgllA2800 (right) for spectral type B2. Meaning of 
symbols/line colors is the same of the previous Figure. 
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Fig. 4. Spectrum of H/3 (left panels) and of MgllA2800 (right panels) for the B1++ source PG 
1201+436. Color coding is the same as for Fig. [2] save for the blue line that traces a blueshifted 
peak visible in both H/? and MgllA2800 fjHand fS.SII . Horizontal scale is rest frame wavelength or 
radial velocity shift from rest frame. Vertical scale is specific fliix in units of 10 



erg s ^ cm ^ 
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Fig. 5. Continuum subtracted spectral ranges of H/3 (top) and MgllA2800 (bottom) 
for RQ (left), CD (middle) and FRII (right) median composites computed in the Bl 

bin. Meaning of symbols is the same of the previous Figures. 
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B1++B1+ B1 B2 Al A2 A3 A4 



Fig. 6. Trends in Pop. A and B as a function of spectral type. The upper panel shows the 
normalized shift the lower S(j). Filled symbols: MgllA2800; open symbols: H/3. 
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-1.5 -1 -0.5 

log W^Edd 



Fig. 7. Trends as a function of Eddington ratio L/L^dd- Upper panel: "dynamical rele- 
vance" indicator = c(|)/FWHM for the MgllA2800doublet broad component. Bottom panel: 
ratio FWHM(MgllA2800)/FWHM(H/3) for the broad components. The smaller circles refer to 
measurements made including the VBC contribution of both lines. 
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3 3.5 4 3 3.5 



log FWHM(H/S) [km s"'] log FWHM(H/S) [km s" 

Fig. 8. FWHM of MgllA2800 vs FWHM H/3 in log scale. Left: f ull profiles, right: o nly BC (of 
single component for MgllA2800). Grey dots are the data points of I Wang et al.l (|2009t ). with the 
line representing their best fit. The large black spots are the median values for the full profiles 
of our median spectra. The dot-dashed lines traces the equality relation, and the thin red line a 
fixed ratio 0.77. 
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